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In  the  previous  paper,  we  presented  a  general  expression  for 
the  Taylor  series  for  any  function  which  can  be  expressed  \s  a 
simple  product  of  a  purely  radial  and  purely  angular  part.  In 
this  paper  we  apply  that  form  of  the  Taylor  series  to  the  Yukawa 
potential  -f-v4a-l-,=„_e  ~  a  v  >  and  to  several  potential  functions 

which  can  be  derived  from  it.  In  particular,  we  examine  the 
expansion  of  the  Yukawa  potential  itself.  In  the  limit  as  the  — 
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exponent  a  vanishes,  the  Yukawa  potential  transforms  into 
the  Coulomb  potential.  We  show  that  the  limiting  process  applied 
to  the  Taylor  scries  for  the  Yukawa  potential  veilds  the  familiar 
form  of  the  Laplace  expansion  of  the  Coulomb  potential. 
Differentiation  of  the  Yukawa  potential  with  respect  to  the 
exponent  a  yields  the  exponential  function.  Hence,  we 
develop  a  Taylor- Laplace  power  scries  representation  of  the 
Morse  potent ia 1  „  1 ntegrat ion  of  the  Yukawa  potential  with 
respect  to  the  enmonent,  in  the  sense  of  the  Laplace  transform, 
yields  functions  \jf  the  form  r'D.  Thus,  Taylor  series  for  the 
Leonard -Jones  and  related  potentials  can  be  constructed, 
finally,  we  consider  a  Laplace- like  functional  expansion  of 
the  Yukawa  potential  followed  by  a  Taylor  scries  develop¬ 
ment  about  the  end-point  of  a  vector  in  the  expansion. 

This  process  illustrates  the  application  of  the  general 
methods  to  a  more  complicated  angular  dependence  than  one 
finds  with  the  use  of  simple  scalar  functions.  An  appendix 
contains  a  Taylor  series  for  one  additional  potential  which 
is  not  directly  related  to  the  Yukawa  potential:  The  Woods - 
Saxon  potential.  The  expansion  for  the  Buckingham  exp-0 
potential  (totally  derivable  from  the  Yukawa  potential!  is 
also  given. 
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Abstract 


In  the  previous  paper,  we  presented  a  general  expression  for  the 
Taylor  series  for  any  function  which  can  he  expressed  as  a  simple 
product  of  a  purely  radial  and  purely  angular  part.  In  this  paper 
w'e  apply  that  form  of  the  Taylor  series  to  the  Yukawa  potential 

.or 

[y(r)  =  e  / r ] ,  and  to  several  potential  functions  which  can  be 
derived  from  it.  In  particular,  we  examine  the  expansion  of  the 
Yukawa  potential  itself.  In  the  limit  as  the  exponent  a  vanishes, 
the  Yukawa  potential  transforms  into  the  Coulomb  potential.  We  show 
that  the  limiting  process  applied  to  the  Taylor  series  for  the  Yukawa 
potential  yields  the  familiar  form  of  the  Laplace  expansion  of  the 
Coulomb  potential.  Differentiation  of  the  Yukawa  potential  with 
respect  to  the  exponent  a  yields  the  exponential  function.  Hence, 
we  develop  a  Taylor-Laplacc  power  series  representation  of  the  Morse 
potential.  Integration  of  the  Yukawa  potential  with  respect  to  the 
exponent,  in  the  sense  of  the  Laplace  transform,  yields  functions  of 
the  form  r  c* .  Thus,  Taylor  series  for  the  Lennard -.Jones  and  related 
potentials  can  be  constructed.  Finally,  we  consider  a  Laplace-like 
functional  expansion  of  the  Yukawa  potential  followed  by  a  Taylor 
series  development  about  the  end-point  of  a  vector  in  the  expansion. 

This  process  illustrates  the  application  of  the  general  methods  to 
a  more  complicated  angular  dependence  than  one  finds  with  the  use  of 
simple  scalar  functions.  An  appendix  contains  a  Taylor  series  for 
one  additional  potential  which  is  not  directly  related  to  the 
Yukawa  potential:  the  Woods-Saxon  potential.  The  expansion  for  the 
Buckingham  exp-fi  potential  (totally  derivable  from  the  Yukawa  potential) 
is  also  given. 


1.  Introduction 


(.liven  a  general  function  G(r)  of  the  vectorial  separation 
between  two  points,  r  =  rj  -  r2 ,  we  can  consider  the  expansion  of  this 
function  in  two  ways.  In  one  form  of  expansion,  we  can  express  the 
original  function  as  a  series  of  products  of  new  functions.  The 
arguments  of  the  individual  functions  depend  only  upon  components  of 
ri  or  r 2  alone.  Thus,  the  dependence  on  rj  and  o  in  the  original 
function  is  separated.  The  other  form  of  expansion  is  to  consider 
the  development  of  the  function  as  a  power  scries  in  terms  of  the 
d i sp 1  a cement ( s )  about  the  vectorial  end-point  (s) .  The  Taylor  scries 
represents  such  a  development. 

In  this  paper  we  develop  Taylor  series  for  several  commonly  used 
potential  energy  functions.  The  functions  all  can  be  derived  from  the 
Yukawa  potential  function: 

y(r]  =  exp(-ar)/r.  (1.1) 

The  particular  functions  we  consider  in  addition  to  the  Yukawa 
potential  are  the  Coulomb  potential  (we  show  this  as  a  limiting  form), 
the  exponential  and  Morse  potentials,  and  the  Leonard -Jones  potential. 
Included  is  the  Buckingham  exp-6  potential  which  is  merely  a  hybrid 
form  of  terms  which  we  consider.  Our  analysis  shows  that  for  parent 
functions  which  can  be  expressed  as  inverse  powers,  the  Taylor  scries 
is  equivalent  to  the  Laplace  expansion  of  the  same  function. 

An  objective  in  seeking  functional  or  power  series  expansions 
of  various  potential  functions  is  to  be  able  to  simplify  the  potential 
energy  functions  for  complicated  distributions  of  molecular  sources. 
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Thus ,  given  a  spatial  distribution  of  sources  for  a  particular 
potential,  e . g .  ,  electrostatic,  exponential,  Morse,  or  lennard- 
Jones,  it  is  frequently  necessary  to  consider  some  form  of  symmetry- 
adapted  expansion.  Uriels  (1980)  has  considered  this  type  of  problem 
for  functional  expansions  of  the  Buckingham  (1938)  cxp-6  potential.  In 
a  number  of  instances,  it  is  necessary  to  consider  symmetry-adapted 
Taylor  series  as  well.  Schmidt,  Pons  and  McKinley  (1980)  considered 
such  an  expansion  in  an  analysis  of  the  vibrations  of  ions  and  atoms 
in  condensed  phases.  The  formulae  we  consider  here  simplify  the 
problem  of  constructing  Taylor  series. 

In  the  next  section  we  list  the  formulae  needed  for  the  subsequent 
discussion.  In  section  3  we  consider  the  Taylor  series  for  the  Yukawa 
potential.  In  section  4  we  consider  the  limiting  transitions  to 
the  Coulomb,  exponential,  and  Lennard- Jones  potentials.  Finally, 
in  section  5  we  consider  the  Laplace  functional  expansion  of  the 
Yukawa  potential  and  the  development  of  a  Taylor  scries  for  one  of 
the  functional  elements.  The  Woods-Saxon  (1954)  and  Buckingham  (1938) 
exp-6  potentials  are  considered  briefly  in  an  appendix. 

1.  The  Ceneral  Taylor  Series 

In  the  preceeding  paper  (McKinley  and  Schmidt,  198_)  we  showed 
that  for  any  function  C(r)  which  is  separable  into  purely  radial  and 
angular  parts  of  the  form 

C,(r)  =  Y^U)I(r)  (2.1) 

where  Y^  (r)  is  the  spherical  harmonic  function,  the  Taylor  series 
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in  terms  of  the  vectorial  displacement  c  is 


G(r+c)  =  (4ii )  3//  2  l  (cu/n!)  [ 


n  =  0 


1, , M ,  £  ,m 


(2L+n 

TTuT)TITTTT 


1/2 

(I.P.OO  |  XO)  (UMm  |  Xp)  I  (r)  . 


(2.2) 


In  this  formula  (L?.Mm|Xp)  is  the  Clebsch-Gordan  coefficient  (Rose, 
1957).  Ay  is  given  by  (Morse  and  Feshhack,  1953) 


A 


nf 


0  for  ?  >  n  and  n  -  C  odd 


(2*+l  )n!  (n- J+l)  !  ! 
7  n  -  7  + 1 )  !  (n  +  ?  +  1  )  !  ! 


for  5.  <  n  and  Z 


n  even 


(2.3) 


and  I  |  (r)  is  defined  by 

on 

dk  kn+  2 f (k) j j  (kr)  (2.4) 

0 

where  in(x)  is  the  spherical  Resscl  function  of  the  first  kind 
(Arfkcn,  1970).  The  quantity  f(k)  is 


Tnh(r)  = 


-A  3 


f  (k) 


.  .  X 
4  n  1 


far 


0 


r2F(r)jA(kr). 


(2.5) 


hquation  (2.2)  and  the  supporting  equations  (2.3)- (2.5)  arc  used  in 
the  following  sections. 

We  shall  have  need  to  use  the  limiting  form  of  cun  (2.2)  for 
a  scalar  function  G:  X  =  p  =  0.  This  quantity  was  derived  in  the 
proceeding  paper  (McKinlcv  and  Schmidt,  198  ) 
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UoU)  =  /'In  ^(c'Vn!  Jh- i)J  +  nAlUl’Jl(r*c)  InJi(r)  .  (  2  .  <>  ) 

!\>(x)  is  the  Legendre  Polynomial. 

3.  Functional  and  Power  Series  for  the  Yukawa  Potential 


We  now  express  the  Yukawa  potential,  eqn  (1.1),  as 

y(r)  =  /47rY0  o  (r)^exp(-ar)  .  (3.1) 


The  Fourier  transform  is 


f(k) 


(4tt)‘5/'2|  dr  r  exp  (-ar)  j  0  (kr) 
0 


(d^/V1 


dr  cxp(-ar)sin(kr) 


(4»rv: 


o 

1 

k 2  +  a : 


(3.2) 


The  radial  factor,  1 n y  C r ) ,  in  the  general  term  of  the  Taylor  series 
i  s 


!n^r) 


dk  k 


n  +  2 


k2 


7  j*(kr) 


(3.31 


The  integral  can  be  evaluated  in  the  upper  half  of  the  complex  plane. 
The  entire  integrand  is  even  (because  n  -  l  is  even).  Thus  the 
range  is  doubled  to  include  the  entire  real  axis  of  k.  We  now 
introduce  the  spherical  liankel  functions  with  the  use  of 
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i  (kr)  =  Re  hj])  (kr)  (  3.1  ] 

n+  ■* 

[as  is  customary,  Re  stands  for  the  'real  part  of'].  The  factor  k 
in  (5.5)  removes  the  pole  of  order  £+1  in  h^^(kr)  at  k  =  0.  Thus,  the 
only  poles  in  the  integrand  are  located  at  the  points  k=Jia.  As 
a  result,  the  integration  yields  modified  spherical  Bessel  functions 
of  the  third  kind  (Arfken,  1970)  which  are  defined  by 

k  .  (a rl  =  i ^  +  “h £  1 ^ ( i a r) .  (5.5) 

Altogether,  we  find 

The  Yukawa  potential  at  the  displaced  point  now  can  be  written  as 

y  ( r + c )  =  l  AnJlan+1P  (r-c)k£(ar).  (5.6) 

n  ,  £ 

It  is  useful  to  compare  this  form  of  the  expansion  with  a  Laplace- 
type  functional  expansion  of  the  same  potential.  [in  section  .5 
we  consider  the  Taylor  expansion  of  part  of  the  following  functional 
form.)  Given  the  Yukawa  potential  y(r),  eqn  (5.1),  we  can  consider 
the  following  relationship: 

>’  (  r  i  -  r  ? )  =  ■— l—  [d  3  k  — exp[-ik*(rl-r2)]  (5.  "I 

(  2  7i  J  3  •'  k  2  +  a  2 


L 
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As  is  shown  in  section  5,  with  the  use  of  two  Rayleigh  expansions 
lor  the  exponential  functions  in  eqn  (0.7),  the  functional  expansion 
i  s 


in  which  i  (x)  is  the  modified  spherical  Ressel  function  of  the 
first  kind.  The  quantities  r>.  ^  represent  the  greater  (lesser)  of 
rj  and  r;.  In  tiiis  expression,  apart  from  powers  of  r  implicitly 
contained  in  the  representat ions  of  the  transcendental  functions 
k  and  i ^  as  infinite  power  series,  there  is  no  explicit  dependence 
on  the  powers  of  ri  or  r2.  The  functional  expansion  is  therefore 
distinctly  different  from  the  Taylor  series.  The  angular  dependence, 
expressed  through  the  Legendre  polynomials,  is  the  only  similarity 
between  the  two  expressions. 


A 

y  ( r  i  -  r  2 )  =  4  7T  l  aix(ar<)kA(ar>)Y  (r 

\,p 


4.  Potentials  Derived  from  the  Yukawa  Potential  and  its  Tavlor  Series 


In  the  Yukawa  potential  y(r),  in  the  limit  as  a  tends  to  :ero 
we  recover  the  Coulomb  potential: 


e~ar 

1 im  — - —  =  1/r.  (4.1) 

a-’-O 


Differentiation  of  the  Yukawa  potential  with  respect  to  the  exponential 
coefficient  a  yields  the  simple  exponential  function: 


e  ~ a  r 
r 


e 


-  ar 


And,  finally,  when  the  Yukawa  potential  is  part  of  the  following 
in  teg ra 1 


TTfW7 


a 


q-2 


(4.3) 


we  recover  functions  of  the  form  r  . 

In  this  section  we  consider  these  three  basic  limiting  cases 
and  the  Taylor  series  which  can  he  derived  from  the  indicated  limiting 
processes . 

The  Coulomb  potential  is  an  appropriate  limit  of  the  Yukawa 
potential,  as  indicated  above.  When  this  limit  is  applied  to  the 
series  representation  (3.61,  wc  obtain  the  following: 


1 

r  +  c 


lim  v(r+c) 
a  -*■  0 


l 

n,f. 


(-cln 

hi 


A  1 im { an+ 1 k  (arl }T\  (c • rl . 
n*a->0  '  * 


(4.41 


From  the  definition  of  the  modified  spherical  Ressel  function 
(Arfken,  1970), 


lim  k  ( a  r ) 
a->0 


lim 

a->-0 


(21)! 


2%!  (ar) £+1 


lim 

a-*0 


( 2  *,  - 1 )  !  ! 


(ar) 


JJT 


(4.5) 


so  that 


1  i  m  [  an+  1  k  (ar)  1  =  p  -ygy  (  2  £  - 1 )  !  !  . 
a->0  '  ’  r  ' 


(4.6) 


From  eqn  (2.3)  which  defines  the  coefficients  A  p ,  wc  find 
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S, ! 

A  f  f  "  XTTTTTT  * 


Ppon  the  substitution  of  eqn  (4.6)  and  (4.7)  into  (3.6),  we  obtain 


I  i;+c ! 


(-c) 

TT+T 


( r • c ) 


(4.8) 


which  is  identical  to  the  Laplace  expansion  of  the  Coulomb  potential. 

We  note  that  if  a  similar  analysis  is  carried  out  with  the  use 
of  eqn  (4.22)  of  the  preccedinp  paper  (McKinley  and  Schmidt,  198  ), 
we  recover  the  Carl  son - Rushbrooke  (1950)  expansion  in  the  limit  as 
the  quantity  a  vanishes. 

The  simple  exponential  function  exp(-ar)  has  been  in  use  for 
a  considerable  time  as  a  representation  of  the  repulsive  interactions 
which  operate  between  atoms,  ions,  and  molecules  (Born  and  Mayer.  1932). 
Combinations  of  the  exponential  repulsion  and  inverse  powers  of 
attraction  arc  commonly  used  in  the  analysis  of  int ermo lecul a r 
interactions  (Hi rsch fel der ,  e_t  ,  1954  ,  Marpcnau  and  Kostner, 

1969).  A  combination  of  exponential  functions  defines  the  Morse 
l 19  29)  potent i a  1 : 


M  ( r )  =  I)  c*(ro-r)[ea(r0-r)  .  2] 


(4.9 ) 


in  which  P  is  an  cnerpy  of  dissociation,  and  r0  is  an  equilibrium 
separation.  Rriels  (1980)  has  considered  the  functional  expansion  of 
the  exponential  component  of  the  Buckinpham  exp-6  potential  in  order 
to  develop  symme t ry - adapt ed  scries. 


The  Taylor  series  for  the  exponential  function  is 


-1  0- 


i  ,  .  n  d  o  x  n  [  -  a  ( r + c )  1 

cxpl-.(ir*c>l  ■  -JIT  I  riel  '  - L 


c  (  -  C  V  ,  n  %  d  ,  11+  1  .  ,  ,  , 

^ .  m  An?,\(r*c)n7Tla  kP(ar)1 

n  .  9 


With  the  use  of  the  mixed  recurrence  relation  (Ahramowitz  and  Stegun, 
19<>5) 


=kj(z)  =  -  zkf  l(z)  -  (  e  + 1 ) k  ?  ( z ) 


(4.11) 


in  eqn  (4.9)  ,  we  get 


exp  [  -  a  ( r  +  c  )  ]  =  £  A  (r  *c)  [ark  f  ]  (ar)  -  (n- 2)k?  (;ir)  |  . 

n ,  ? 

(4.12) 


The  Morse  potential  expands  as 


H<ro  -  i)  I  ■ 4^4-  An( 


Pp  (r-c)  {2ne“ar°(  2ar  k?  ^  (2ar) 


(4.131 


(n- 9.)k?  (2ar)  1  -  2  car°[ar  k?1(ar)  -  (n  -  Z)  (ar )  1  } . 


finally,  we  consider  the  Tavlor  series  for  | r+c |  .  Substituting 

the  Taylor  series  for  ttie  Yukawa  potential  (3.6)  into  (4.3),  we 
obta  i  n 


r+c i ^  ' ^nr  j; 


*  go 

rTTTTT  J  d: 


1  *q'2 


n  f. - nT~  1  y  kx-.c  j  j 


Pn  (r*c)  f  da  ;/1  +  n  1  k 0  (ar ) 


0 


(4.14) 
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where  we  have  exchanged  the  order  of  summation  and  integration.  Ke 
now  can  use  the  mixed  exponent ial -polvnomi al  form  of  the  modified 
spherical  Bessel  function  (Abramowitz  and  Stegun,  1905) 


(  £  +  s  )  ! 

(2s')!!  ( l  -  s )  !  x 


(4.15) 


to  write 


da  ac,  +  n  ^k?(ar) 


1  9 

7  l 

r  s  =  0 


(g+s)! 
(2s)  !  !  ( \  -  s  ) 


1  ^ 

>J. 


,  q  +  n-s  -2  -ar 
da  a  1  e 


(4. 1(0 


The  integral  is  uncomplicated  and  straightforward  to  evaluate  as 
long  as  q^2 .  (Note,  from  eqn  (4.14)  the  expansion  vanishes  for 
(q-2)<0  because  the  factorial  becomes  infinite.]  Thus,  the  expansion 


i  s 


l  AnPP/r<C) 


.  i  q  q  t  . '  n  £  9 
r  +  c I  r  1  n  ,  £  ,  s 


(J  +  s)  !  (c 

1  +  n- s  -  2 )  ! 

n  !  (  £  -  s  )  !  ( i 

-  s  )!  !  ( q -  2  )T 

s  f o rm  is 

equi va 1 ent 

obtained  by  Briels  (1980)  who  used  a  functional  expansion 


5.  Laplace  and  Laplacc-Taylor  expansions  of  the  Yukawa  potential 

The  main  problem  we  consider  in  this  section  concerns  the 
evaluation  of  the  Taylor  series  expansion  of  the  function  y(rj-r2) 
about  the  point  ri.  Such  an  expansion  is  required,  for  example, 
when  one  needs  to  consider  a  source  at  r2  and  the  displacement 
of  a  particle  about  the  point  ri;  both  points  r,  and  r2  share  a 
common  origin. 
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Bcfore  wo  consider  the  displacement  about  rj  by  c  in  detail, 
we  first  demonstrate  eqn  (3.8),  the  Laplace  functional  expansion 
of  the  Yukawa  potential. 

The  louricr  transform  is  given  by  eqn  (3.2)  together  with  Y00(k): 


!■' .  t.l  y  (  r)  } 


(5.1) 


where  here  q  is  the  strength  of  the  source.  The  functional  form  of 
)’  ( r  i  -  r  2 )  is 


y  (  r  i  -  r  2 )  =  -  1  fd3k  -  -  T- L-  exp  f  -  ik*  (r  i  -  r2 )  ]  . 
(  2 tt  )  3  J  k 2  +  a  2 

The  use  of  two  Rayleigh  expansions 

exp  (ik*  r)  =  4tt  l  i  X  Y*  ( r)  Yx  (k)  j  A  (kr) 

A.u  "|J 


(5.2) 


(5.3) 


in  eqn  (5.2)  yields 

y ( r ! - r 2 )  =  8q  £  vAp (r i ) YAp (r? )J  dk  iA(kri)jA(kr2)  (5.4) 

A  ,  u  J  o  k  +  a 

upon  integrating  over  the  angles. 

The  radial  integral  is 


1 (ti ,r2) 


dk - -  j  v  (kri) j ,  (kr 2 ) . 

o  k2+a2  A  A 


(5.5) 


The  evaluation  of  this  integral  proceeds  along  lines  similar  to  the 
evaluation  of  eqn  (3.3).  We  write 


13- 


I  =  ike 


Jk  7777  h»1)(kr>'j»‘kr<) 

,»  k^+a' 


(  5  .  (>  ) 


in  which  r>  is  the  greater  of  ri,  r?.  The  (.A+l)-order  pole  in  Ik* 
at  k=0  is  removed  by  the  (A  +  2) -order  zero  in  k2j.(kr,).  lienee, 
closing  the  contour  above,  the  only  pole  is  at  k=+ia.  The  integral 
l  is 


1  =  7aix(ar<)kyCar>)  • 


5.7) 


Thus,  the  Laplace  expansion  of  the  Yukawa  potential  is  that  given  by 
eqn  (.3.8)  : 


y  ( r  i  -  r  2 )  =  4  ti  q  £y  Au( r  i  )  Y  A  (r  2 )  a  i  A  (ar  <)  k  (a  r  > )  . 
Ap  M  M 


(3.8) 


We  now  consider  the  expansion  about  the  point  ri 


00  1 

y  (r  i+c-r2)  =  q  £  i,  (c  •  V,)  ny  (r  x  -  r  2 ) 


n  =  0 


=  4iu'  l  iY'YIM(r2)(~’^n fYxMtri)aix(ar<)kx(ar>)  I  • 

nAp 

(5.8) 

In  order  to  evaluate  this  expansion,  we  need  to  consider  the  'partial 
potential ' 


=  YXu(ri)aiX(ar<)kX(ar>) 


1 2n) 


d 3 k  fA(k)YA^j(k)cxp(-ik*r1) 


(5.0) 
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whc  re 


MM  = 


,  -Af 
4  ii  i 


dri  rt4'x(ri)3x(kri) 


-  4  it  i 


.  A 


f 1  - 


dr  !  rta i , tar j  Jk.  (ar2 )  j  .  (.kr j  ) 


15.10) 


+  d  r  2  r  j  ai.(ar2)k,(ari)j,(kr. 


We  show  in  appendix  2  that 
.  -  A 

1  rAlM  =  ;j  A(kr2J. 

A  k2  +  a2  A 

In  the  standard  manner,  wo  now  find 

il:  u  '  •  ’ " ;  Or‘’  •  0»)3/2  £  E  (-i)"‘\tYt.(Ovult;,) 

fcmLM 

{  ’i  +  I  1  1/2 

x  [  i:\+T)Tl 4+0  J  U^OO  I  AO)  (LtMro I  A,.)  InLX  ( r , ) 


(5.11) 


(5.12) 


where  specifically. 


Wro  = 


12.)  * 


,  n+2  4ni A  ...  ^ 

dk  k  -  j.  (kr2)  1,  (kri) 

k2  +a 2  A  L 


.  A 


2  n  2 


,n+2 

dk  — - j  .  (kr2)  j  ,  (kr  ]  )  . 

,  k 2  +  a 2  A  L 


(5.1 3) 


from  the  parity  of  the 
Anf,  n-£  =  even.  Thus, 
doubled.  Write 


C, l ehsch- (Iordan  coetfic  ient 
h+n+A=cven.  The  integrand 


4+L+\=oven.  from 
is  even  and  can  he 
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nLA 


Ui) 


dk 

-  oo 


kn+: 

- jA(kr2)j.  C k r i )  . 

k2  +a2  A  L 


(5.14J 


We  consider  two  cases  separately. 

Case  A:  rj<r2.  Let  j^(kr2)  =  Re  hj^-*(kr2)  and  close  the  contour 
above.  As  k-+0,  the  entire  integrand  behaves  as  +  The 

The  triangle  inequality  for  the  Clebsch- Gordan  coefficient  ensures 
that  £+L-A>_0.  We  know  that  n-Jl^O.  Thus,  n+L-A>0,  and  there  is  no 
pole  at  k=0.  The  pole  at  k=ia  contributes  through  its  residue  to 
the  integral: 

-n  +  L  . 

InUlri)  =  —  a”  1kx(ar2)ix(ar1)  (5.15) 

4  7T 


and 


y(ri+c-r2)  = 

-  ~  r  i  <  r  2 


2L+ 1 

(Ta+1)  ( 2 £  +  1T 


(47r)3/2q  l  (cn/n ! ) A  . Y*  (r 2 ) Y ?  (c) Y 


nAy  £mLM 


n£  Ay 


Slm^'  LM 


1/2 


(LfcOO  I  AO)  (StLmMl  Ay)aI1  +  1k,  (ar2)  i  (ar,). 

A  L 


(ri) 

(5.16) 


In  the  limit  as  rj+0  in  eqn  (5.16),  we  recover  eqn  (3.6). 

Case  B:  ri>r2.  Proceed  in  a  manner  which  is  similar  to  case  A. 
Let  j  j  C  k  r  i )  =  Re  h^^(kri).  Again,  wc  can  show  that  n  +  A-L^>0  so  that 
there  is  no  pole  at  k=0.  The  residue  at  k=ia  now  yields 


nLA 


(r  i) 


. n  +  2  A  - 1 
4  7 r 


an+liA(ar2)kA (ar,) 


(5.17) 


and 
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>'(r,*c-r,)ri>r2  -  (4»)3/2q  J  (cn/n  0  A„ ,  Y?„  ( r,  )  V,m  (O  Y,  „  (  r ,  ) 


nApfcmLM 


nf  An1'  2  "  '  1,M 


?i+l  \  1  /  2  , 

x  1  nT^TTTTrnrj  (LAOO|  ilO)  (Lil  Mm|  Ap)an  iA(ar2)kL(ar,)  .  (5.18) 


As  an  example,  we  now  apply  these  results  to  the  problem  of  a 
ring  of  source,  such  as  a  charged  ring  or  a  ring  of  source  for  the 
Morse  potential. 

The  distribution  of  source  is  such  that  an  element  of  source  dq 
is  in  the  vicinity  of  r2.  Thus,  we  write 


y(i'i)  =  dy  i  -  r  2  ) 


=  4nfdq  l  Yx  (r^Y*  (r2)  ai^ar^k^ar^ 

J  A,p 


(5.19) 


for  the  case  of  a  ring,  the  magnitude  of  r2  is  constant.  Hence, 


dq  =  d<j); 


(S.20) 


’Vr’> 


q  Yxo^r2^6M0‘ 


(5.21) 


We  find  for  y(ri) 


y(ri)  =  qa^>A(ar<)kA(ar>)  (2A  +  1)Pa(cos0j)Pa(cos62) 
A 


(5.22 
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Wc  now  use  the  results  obtained  above  to  investigate  dis¬ 
placements  by  c  about  r i.  For  ri<r2  we  have 


'Oi*c)ri<  -  (4.)3/2q  I  ii(cn/n!)An/>0tr,)YLH(r,)Y  (O 


nAp  2.mI,M 


1/2 


T7x^nT~HTr)  ' '  "o^00|A0)(L£Mm|X0)an  +  1i1  (ar!)kx(ar2)  (5.25) 


and  for  i"i>r2 


(ri  +  c)  =  (4tt)3//2(i  l  — C  ■ 


2L  +  1 

T2T+ 1TC2  £+T) 


nAy£mhM  n! 
1/2 


AntY»0<rOV,,M(r,)Ylnl(c) 


(LtOO|i())(L«Mm|XO)a"41i1(ar2)kI  far,).  (5.24) 


In  the  limit  as  a  tends  to  zero  in  these  expressions,  we 
recover  the  forms  associated  with  the  Coulomb  potential: 


1  ini  y(ri) 
a->0 


q£(2A+l)P,  (cos0i)R.(cos02)lim  ai.  (ar  ) k.  (ar  ) 


a-*0 


q^Px(cos0i)Px(cose2) 

A  r. 


(5.25) 


from  eqn  (5.22).  This  expression  usually  is  obtained  as  an  example 
in  classical  electrostatics  by  other  means  (cf . ,  Jackson,  19b2,  p.b4). 
The  Coulomb  limit  for  the  expansions  (5.23)  and  (5.24)  yields 


1  im 
a-*0 


yfr»+£-r^ri<r2 


( 4 7T )  3,/2q 


l 

?  km 


T!Yk-H),n  ^r2^YL 


-m 


(r, 


>Ylmfc> 
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1/2 

.TTFTrrTrWrr)  at»»u^»)at-Miiuto).^^||j- 


A 

<_T+T+T 

r2 


(5.26) 


and 


3/2  r 

lim  y(ri+c-r2)  =  (4tt)  q  £ 

a+0  ~  ~  ~  r ! >r2  lLm 


xY«n,<c> 


2L  +  1 

U2H+1)  (2i-2L+l) 


1/2 


(LJ.00  |  £  -  L  0 )  (L$,-mm  f  S.-LO) 


!>-L 

( 2  L - 1 ) ! !  r2 
(21-2L+1) ! !  L+ 1  ' 


(5.27) 


Although  we  have  examined  the  case  of  the  ring  as  an  example, 
it  is  far  from  an  idle  one.  There  are  numbers  of  systems  of  interest 
in  physics,  chemistry,  and  biology  in  which  essentially  one  atom  or 
ion  sits  in  close  proximity  to  a  configuration  of  atoms  which  are 
bound  together  as  pentagons,  hexagons,  and  higher  regular  polygons. 

Some  examples  are  the  single  crystal  surfaces  of  pure,  clean  metals 
and  other  solids,  annular,  charged  ring  molecules  in  chemistry, 
and  surface  aggregates  of  phosphate  ions  in  biological  membranes. 

It  is  easy  to  show  by  computer  simulation  (Schmidt,  unpublished) 
that  the  difference  between  the  hexagon  and  a  continuous  ring 
is  a  small  one.  Thus,  the  ring  of  charge  or  matter  as  the  source 
density  for  the  Morse  potential,  for  example,  serves  a  useful  function. 
The  analyses  of  vibrations  and  stabilities  in  such  systems  are 
facilitated  by  the  analyses  given  above. 


m 
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Appendix  I.  The  Woods-Saxon  and  Buckingham  exp-6  potentials 

We  have  mentioned  the  Buckingham  (1938)  potential  in  the  text 
of  this  paper.  In  this  appendix  we  show  a  general  expression  for  the 
Taylor  series  expansion  of  this  potential.  The  result  compliments 
Bricls  (1980)  functional  expansion.  first,  however,  we  consider 
a  different  kind  of  potential  energy  function,  one  which  depends 
upon  the  exponential  function,  but  is  not  directly  related  to  the 
Yukawa  potential.  The  function  is  the  Woods-Saxon  (1954)  potential 
which  has  been  used  frequently  in  the  analyses  of  nuclear  models. 

The  Woods-Saxon  potential  has  a  simple  form: 

^WS  -  1  +'  exp  [  (r-r0)/p]  ’  U-l) 

This  function  does  not  easily  admit  a  Fourier  transform.  Therefore, 
the  direct  application  of  the  integral  form  of  the  Taylor  series 
is  inappropriate.  We  resort  to  a  differential  form  (McKinley  and 
Schmidt,  198  )  in  order  to  get  useful  results. 

In  particular,  in  the  previous  paper  we  showed  that  for  a 
scalar  function  of  the  form 

G(r)  =  Y„  o  (r)  [  /4ttF  ( r  )] 

=  F(r)  (I  .2) 


an  arbitrary  term  in  the  Taylor  series  is 


1  fr.n1nf:frl  _  r  a  r  (-l)qU+q)i  -qfn-q 

n!Ct:  V  ST  £AnsT  3qi0  r  (  r 
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+ 


(d/dr)n"q_1F(r) . 


(1.3) 


For  the  case  of  the  Woods-Saxon  potential,  it  is  reasonably 
easy  to  obtain  a  closed,  polynomial  representation  for  the  differentia¬ 
tions  indicated  in  eqn  (1.3). 

The  differentiation  of  the  Woods-Saxon  potential  (1.1)  to 
arbitrary  order  is  carried  out  as  follows.  Let 


x  =  exp [ (r-r0)/p]  . 


(1.4) 


Then 

(d/dr)n  =  P'n(x^]n  (1.5) 

and  by  means  of  mathematical  induction  we  find 


11  =  l  cV(d/dx)s  (1.6) 

s  =  l 

where  the  coefficients  are  determined  by  means  of  the  following 
initial  and  end  conditions 


C?  =  Cn  =  1 
1  n 

Cn  =  0  for  all  s  >  n 
s 

and  the  recursion  relation 


(1.7) 


(1.8) 
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The  differentiation  of  the  Woods-Saxon  potential  now  leads  to  the 
general  term  for  the  Taylor  expansion 


1 

n ! 


,  U-v) 


%  =  -  Voc^  yA  p  rC.c)  y  fD/r.q 

45  ws  _,.n  rnr  £lr  c}  An  RIW!  tP/  J 


n !  p 


n  -  q  - 1 


q  =  0 


U-q)(p/r)Y'  e  (s)  +  "jV'*1  e(s) 


s  =  1 


s  =  l 


(1.9) 


in  which  e(s)  is 

els)  -  (-l)s  - ..£*£[  gllll.  9 I/-£i —  .  (i.io) 

{1+exp [ (r-r„)/p] }b 

The  Buckingham  potential  is  simply  a  combination  of  the  exponential 
and  van  der  Waals  r  b  potentials: 


4>R(.r)  =  ae‘br  -  c/r6. 


(1  .11) 


The  Taylor  expansion  of  this  function  is  just  the  combinations  of 
eqn  (4.12)  and  (4.17).  We  find  specifically 


<t>B(r  +  r  ’) 


1 

U 


j  —  — '  1  A  P  (r*r') 
L0  n!  rtn£tVr  r  ; 
n ,  £ 


abn[brk£ _  j (br)  -  (n- £) k£ (br)  | 


c  v  (£  +  s)  ! (n+4  - s )  !  ] 
rn+6  s^0  (£-s)  TT2JTTT  )  * 


(1  .12) 
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Appendix  II.  The  evaluation  of  the  integral  (5.10) 

We  evaluate  the  Fourier  radial  transform,  eqn  (5.10), 
as  follows. 

Write  the  spherical  Bessel  equation  as 

a7[r23?jx(kr)l  +  (kr)2jx(kr)-  X(X+l)j x(kr)  =  0.  (II. 1) 

We  write  the  modified  spherical  Bessel  equation  for  any  of  its  solutions 
as 


3?tr2a7Mar)]  '  C“r)2SA(ar)  -  X(X  +  lKx(ar)  =  0.  (II. 2) 

Multiply  the  first  equation  by  -^(ar)  and  the  second  by  jA(kr)  and 
add: 


^“^aT^crrj  X^kr^  +  J A(kr) [r2jpCx(ar) ]  -CA(ar)  (krftx(kr) 


-  jA(kr)  (ar) 2tA(ar)  =  0 


(II. 3) 


The  first  two  terms  are  integrated  by  parts.  The  third  and  fourth 
terms  are  transposed.  The  result  is 


r2C,(ar)j  (kr)  =  — - -  g- 

X  X  k2+a2  dr 


d  fjx(kr)r2^x(ar)  '  r2cA(ar^jx(kr) 


(11.4) 


When  this  result  is  substituted  into  eqn  (5.10),  we  find 


Vk)  =  rf77(akxCar2)[UCkri)r2l3F7ix(ari)  •  r2ix(ar^HTTjx(kri-l] 


k2+a 


+  aix(ar2)  [jx(kr1)rf^-kx(ar1)  -  r?  kx  (ar  i  )^4-j  x  (kr  i )  ]  0 

1  1  r 


A  'A 

~~~  a2r2jx(kr2)  [kx  (ar2)  ij[  (ar2)  -  ix(ar2)k^  (ar2)].  (H.5) 

k  ^  ct 


The  quantity  in  the  brackets  is  a  simple  Wronskian  which  has  the  value 


l/(a2r2)  (Arfken,  1970).  Thus, 


fX^  = 

k  +a2 


(5.11) 
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